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We extend the Exchange Fluctuation Theorem for energy exchange between thermal quantum 
systems beyond the assumption of molecular chaos, and describe the non-equilibrium exchange 
dynamics of correlated quantum states. The relation quantifies how the tendency for systems to 
equilibrate is modified in high-correlation environments. Our results elucidate the role of mea- 
surement disturbance for such scenarios. We show a simple application by finding a semi-classical 
maximum work theorem in the presence of correlations. 
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Fluctuation theorems describe non-equilibrium trans- 
formations of thermodynamic systems and constitute a 
refinement of the second law of thermodynamics, the 
most well-known incarnations being the work fluctuation 
theorems due to Jarzynski and Crooks [H-Q- However, 
in addition to focussing on the extraction of mechanical 
work from a single quantum system, an equally funda- 
mental topic is the thermodynamic tendency of multi- 
partite systems to equilibrate. The canonical example 
of this is heat exchange between two thermal systems at 
different temperatures and leads us instead to fluctuation 
theorems for heat that provide a quantitative description 
of the fluctuations in energy exchange between two hot 
bodies. 

The thermodynamic arrow [H, Q is one particular man- 
ifestation of the second law of thermodynamics, and in 
its canonical form states that on average, heat will flow 
from a hotter body to a colder one. Specifically, given 
two thermal states pA and at temperatures Ta and 
Tb with respect to Hamiltonians Ha and Hb, and an 
energy-conserving unitary evolution of the joint state, 
PA® Pb U PA ® PbU\ we define heat flow into A as 
Qa = ^t:[Hap'a\ ~ tr[i?APA] where is the final reduced 
state for A, and we assume that the free Hamiltonians for 
A and B do not change. The fact that Gibbsian states 
minimize the free energy yields the inequality 



Qa{t^ 7^) 

J^A J^B 



>0, 



(1) 



and so if Ta < Tb we have that the Qa is strictly posi- 
tive, and on average energy is transferred from the hotter 
body to the colder one. This is the standard thermody- 
namic arrow for heat flow. 

However, a sharper expression of the directionality for 
heat flow exists in the recent Exchange Fluctuation Theo- 
rem (XFT) due to Jarzynski and Wojcik which states 
that for two systems A and B, initially at temperatures 
Ta and Tb, the probability P{q) of a sharp exchange of 
energy q from B to A obeys the relation 



P(q)/P(-<z)=cxp[A/3g], 



(2) 



where A/3 = (fcr^)^^ — [UTb)^^ ■ This relation quantifies 
the relative likelihood of a fixed exchange process and 
its time reversed twin, and shows that heat fiow from 
a colder to a hotter object is exponentially suppressed. 
A simple application of Jensen's inequality to ([2]) leads 
to an averaged inequality {q){l/TA — ^/Tb) > 0. This 
seems to suggest that ([T|) automatically follows from ^ , 
however it must be emphasized that while Qa equals 
(q) for classical states, this need not be true for more 
general quantum mechanical states. In general, a sharp ( 
rank-1 projective) energy measurement will produce non- 
classical disturbances in the quantum state of a system, 
and even though the expression (g)(l/TA— 1/Tb) > still 
provides a thermodynamic directionality, it is no longer 
identical to ([T]), which does not involve any measurement 
disturbance. 

The scope of the XFT is impressive, being valid for ar- 
bitrary unitary interactions between A and B that con- 
serve energy, and rests crucially on the two assumptions 
of (I) initial Gibbsian states, and (H) the assumed time- 
reversal invariance of the underlying dynamics. However 
the directionality for the thermodynamic processes relies 
on the third assumption (HI) that the systems involved 
arc initially uncorrelated; namely Boltzmann's assump- 
tion of 'molecular chaos' [?[. As such, both the Clau- 
sius relation ^ and the Jarzynski- Wojcik Fluctuation 
Theorem ([2|) arc limited in application, and will fail to 
hold within the domain of high-correlation environments. 
Indeed, with the extremal case of a globally pure, mul- 
tipartite quantum state with thermal subsystems there 
should exist no directional constraint whatsoever, and for 
such situations no equality such as ([2]) should hold. The 
reason for this is that pure quantum states are states of 
maximal knowledge and may be reversibly intcrconverted 
through the appropriate unitary transformations. Such 
pure state thermality is of importance and turns out to 
be the typical scenario with respect to the Haar measure. 
Specifically, for a randomly chosen multipartite state j^") 
of a system with fixed energy, the reduced state for a 
small subsystem is exponentially likely to be Gibbsian 
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[10|, with the thcrmahty arising due to quantum entan- 
glement. 

Quantum correlations can be far stronger than their 
classical counterparts, and in addition to the Gibb- 
sian typicality in pure states, entanglement theory has 
other deep connections with thermodynamics often 
through their parallel formulations as resource theories 
[T^ . Beyond the foundational interest of studying the 
dissolution of the thermodynamic arrow due to strong 
correlations, there is also rapid experimental progress in 
the precise manipulation of small quantum systems de- 
signed to function as engines at nanoscales, and as such, 
it is also of practical importance to determine the fun- 
damental limitations and behaviour of heat exchange in 
such quantum systems. 

The purpose of this letter is to remove the third fun- 
damental assumption (III) of molecular chaos, and ex- 
tend the existing XFT into high-correlation environ- 
ments, in which initially correlated quantum systems are 
allowed to evolve under non-equilbrium dynamics and 
exchange heat. In doing so we identify the appropriate 
information-theoretic measure for the effect of correla- 
tions on sharp energy exchanges and describe how it may 
contribute in a maximum work theorem scenario. We 
also demonstrate the strong constraints that measure- 
ment disturbance imposes on the formulation of a more 
general fluctuation theorem that allows POVM measure- 
ments on the systems. 

Details of the generalized setting. We wish to establish 
an exchange fluctuation theorem for a bipartite quan- 
tum system, whose reduced states pA and ps are ini- 
tially thermal, but where we drop the assumption of the 
initial factorization of the joint state pab and allow gen- 
uine quantum mechanical coherence and entanglement to 
either evolve or be initially present. 

In defining heat exchange in Q, the isolated bipar- 
tite system is assumed to undergo a three step process, 
in which an initial energy measurement A^i is first per- 
formed on the two subsystems which are then allowed to 
subsequently interact and evolve under a unitary [/, un- 
til a final energy measurement A^2 is performed. Thus, 
the bipartite quantum state pab undergoes the following 
sequence of quantum operations: pab -MiIpab] 
UoMiIpab] M2 o U o Mi[pab]- 

As mentioned, the central assumptions of the XFT 
are time-reversal symmetry of the underlying dynam- 
ics and the initial thermality of the individual subsys- 
tems. Specifically, we assume that Q'^HaQ = Ha 
and &^Hb& = Hb, where & is the anti- unitary time- 
reversal operator. In what follows we use joint en- 
ergy eigenstates \(j),x), such that x) = E^\(l),x} 
and Hb\(I)tX) = for which we deduce that 
ifA0|0, x) ^ E^Q\(l),x) (with a similar expression for 
B). We take {\(f>)} and {|x)} to be complete orthonor- 
mal bases for A and B so that Ha = J2,p-^<p\'f'){'t'\ ^'^'^ 
Hb = J2x^x\x){x\j and since 9 is a symmetry of clas- 
sical states we assume that 6(|0) (S> \x)) is always in 
the basis set {|0) (E) \x)} for any (j) and x- The ther- 



mal marginal states of pab are then given by pA = 

^I'E^e-^^^^I^X'/'l and pB = e-''"''Hx> (xl, 

where Za and Zb are the usual partition functions for A 
and B. 

The measurements Aii and A^2 used to determine the 
energies of the subsystems can in general be POVM mea- 
surements, however we restrict both the initial and final 
measurement to be rank-1 projective measurements onto 
the local energy eigenbases of A and B; a point we return 
to later. Given this, it is useful to introduce the notion 
of a history for the composite quantum system AB. A 
history is denoted 

7 - [pab; \4>) ® \x) A 10') ^ Ix')) , (3) 

where pab is the initial quantum state that is first pro- 
jected into the energy eigenstate |0) (8) |x) under Mi and 
then evolves unitarily to C^(|</>) (X" |x))i which is then mea- 
sured and projected into the energy eigenstate |0') (8) \x') 
under A^2- 

We denote by F the full set of all histories {7} com- 
prised of first beginning in the state pab, measuring out 
some energy eigenstate, evolving under some U and then 
measuring out some final energy eigenstate. 

The thermodynamic condition of energy conserva- 
tion we use is simply that tT[pABiHA + Hb)] = 
tippABU^HA + Hb)]- We note that U involves the 
interaction Hamiltonian that is assumed to be smoothly 
switched on and off, but the energies of the subsystems 
are always measured with respect to their appropriate 
free Hamiltonian. 

While the exchange fluctuation theorem is a refinement 
on the thermodynamic arrow for heat flow, in that it 
deals with subensembles of postselected outcomes with 
sharp energy transfers q, not all of the uncertainty is at- 
tributable to the statistical mixture of different energy 
states. Pure quantum mechanical states allow the possi- 
bility of intrinsic quantum fluctuations, and so while we 
might step beyond classical statistical fluctuations by fo- 
cussing on individual pure state outcomes in the XFT, 
we might also allow the possibility of quantum coher- 
ence evolving under the unitary dynamics and generat- 
ing new indeterminacy. For example, with respect to 
the statistics of energy measurements, the pure quantum 
state lip) cx J2k ^ e~'^-^'» |i?fc) is indistinguishable from a 
thermal mixed state. 

The energy measurement of such superpositions will 
display quantum fluctuations, some of which may in- 
crease the total energy of A_B, some decrease it, but on 
average no net energy should be gained from the fluctu- 
ations. It is then simple to allow histories with positive 
energy fluctuations that increase the total energy of A_B, 
or negative fluctuations that decrease the energy of AB^ 
and also histories with no fluctuations at all. As such, a 
useful and physically intuitive division of the set of histo- 
ries is into sets of histories with similar energy transfor- 
mations. In particular we write F = Ug,A£F((7, Ae), where 
r((7, Ae) is the set of 7 (of the form ^) with fluctuations 
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of the total energy of AB 

Ae^{(t>\x'\{HA + HBM',x')-{^,x\iHA + HBM,x), 
and with an energy transfer q into A defined a,s q — 

An XFT for correlated quantum states. Within the 
setting just described we can now formulate an XFT that 
drops assumption (III) of molecular chaos. 

For an initial state pab the occurrence of a single his- 
tory 7 e r{q, Ae), has probability 



Prob[7] = {cl),x\pAB\q^,x)m',x'),U\q^,x)W 



(4) 



where U = e~'-"* and H = Ha + Hb + i^int is the total 
Hamiltonian, including the interaction i/jnt between A 
and B, which is switched on at the initial time, and for 
clarity we write the Hilbcrt space inner product as (■,■)■ 
From time-reversal invariance, Q'^HQ = H, and the 
anti-unitarity of 6 it follows that U = Q^U^Q, and so 
we have that 

m',x'),u\<j>,x))\' = |(e|0,x>,c/e|0',x'»P- (5) 

Thus, time-reversal symmetry alone implies the probabil- 
ity to go from the initial state \(j),x) to \(j>',x') is always 
equal to the probability to go from Q\(j>',x') to Q\(j),x) 
under the same unitary interaction. In what follows, we 
use a star to denote time-reversed objects, for example 

|<^',x')*:=e|f,x')- 

To quantify correlations in the quantum state as they 
relate to the XFT we define, for any joint local POVMs 
{Mi} on A and {Nj} on B, the quantity I{pAB] Mi, Nj) 
via the expression 

Tf M 1 f tY[M,®N^pAB] \ 

IipAB;M.,N,) := ^-[^,,[M.pA]tr[N,pB]))- 

For a sharp enery measurement, we simply have Aii = 
{M^<E)N^} where {M^ = and {Af^ = \x)ix\} are 

the rank-1 projectors in the energy eigenbases jl8| , and 
when PAB is a correlated quantum state having thermal 
marginals, we find that Mi maps pab into a classically 
correlated state, diagonal in the energy eigenbasis. More- 
over, it is readily seen that tr b{A) [Mi[pab]] = Pa(b), and 
so the state A^i[pas] has the same thermal marginals as 

PAB M- 

For this particular initial measurement, the probability 
p(|0, x)) of projecting into the state |0, x) under Aii can 
be written as 



P(I'/',X» 



,-l3AE^-l}BE^-\oti{ZAZB) + 'I[pAB\M^,N^\ 



(J) 



while a comparison with the probability of obtaining the 
outcome implies that 



(8) 



where A/3 = Pa — 13 b, and crucially AI(7) = 
I[pab;M^>^,N^>J - I[pab; M^, N^] is the appropriate 



correlation measure, dependent only on the initial state 
Pab and the initial measurement Mi. 
Combining ([S]) with ([5]) we find that 



Prob[7] 



„A/}g+/SBA£-AI(7) 



Prob[7*] 

where 7* is the time-reversed twin of 7 given by 



7* = (pab; 10', x')^ 



(9) 



(10) 



In particular, the history 7 involves a quantity of en- 
ergy q being transferred into A and a net increase of 
total energy Ae, while 7* involves the opposite changes, 
and so T{q,Ae)* = T{-q,~Ae). We also note that @ 
is independent of the specific form of the dynamics (be- 
yond time-reversal invariance), and depends solely on the 
properties of the initial quantum state pab- 

One can now compare the ratio of probabilities of the 
set r(g, Ae) and its time-reversed twin set r(g, Ae)* = 
T{—q, — Ae). The probability of the former is given by 



Prob[r[q, Ae]] 



ier{q.Ae) 



A^,+feAe^AI(7)p^0b[7*].(ll) 



While the term g^'^^+'^BAe j^ay be factored out of the 
sum, the correlation term cannot as it will generally 
vary over the set T{q,Ae). Instead we necessarily ob- 
tain bounds for the ratio of the probabilities. To fix the 
lower and upper bounds, we respectively define AI/ = 
max^gr(q,Ae)[AX(7)] and AIu = min^gr(q,Ae)[AX(7)], 
and immediately deduce that 

Al3q+I3BA€-Ali < Prob[q, Ae] A/Jg+feAe-AI„ 

- Prob[-q, -Ae] " ^ ^ ' 

Equation ([T2|l is a generalization of ([2]) and is a con- 
straint on the relative likelihood of a forward transition to 
a backward transition given an initially correlated quan- 
tum state with thermal subsystems. As can be seen, by 
increasing the degree of correlations present in the initial 
state one gradually weakens the constraint on the ther- 
modynamic arrow, as expected. Moreover, it is not pos- 
sible to tighten these bounds without making additional 
assumptions as to the particular form of the dynamics. 

Beyond the relative likelihood of the forward and re- 
verse processes, one can take and sum over 7 G P, to 
obtain the non-equilibrium equality for an initially cor- 
related state 



-AI3q-l3BAe+AI\ 



1 



(13) 



and then using Jensen's inequality we have that 

A/3(g)-f/?B(Ae)-(AZ) >0. (14) 

Here, {AT) represents the difference in the classical mu- 
tual information of measurement outcomes between the 
initial and final states. 
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Given the assumption that (Ae) = we have A(3{q) — 
(AX) > 0, which reduces to the Clausius relation 
Q{1/Ta - l/Ts) > for Q = {q) and the assumption of 
molecular chaos. More importantly it displays the ener- 
getic value of correlations in providing a modified lower 
bound of Q(1/Ta - l/Ts) > (AI) with the function 
^{pab't Mi, Nj) as the appropriate sharp-outcome mea- 
sure for the initial bipartite quantum state. This must be 
compared with averaged results obtained previously d, 
in which (AX) is replaced with the change in the quan- 
tum mutual information of the state pab- The origin of 
the difference is that the XFT demands sharp energies at 
the initial and final stages, as opposed to bluntly looking 
at expectation values of energy for pure quantum states. 

While it is natural to impose energy conservation, ei- 
ther at the level of commuting Hamiltonians or expec- 
tation values, the XFT given by equation (|T2l) makes 
predictions for a particular type of state with local tem- 
peratures and global correlations, and provides only the 
relative likelihood of seeing one forward thermodynamic 
process compared to its reverse - not whether it occurs 
at all. As such, the specific interaction Hamiltonian that 
is used only serves to predict the absolute likelihood of 
these different processes. 

Beyond sharp energy measurements. As mentioned, 
the sharp energy measurements A4i used are quite de- 
structive of coherence, and so one might wonder whether 
an XFT can be obtained for more gentle POVMs. In 
other words, can wc perform the time- reverse pairing 
trick using mixed quantum states? 

Given a preparation of some ctab by the initial mea- 
surement Ail, we wish to do the pairing trick with the 
state aAB and a time-reversed twin. If we drop the as- 
sumption that A^2 is a sharp energy measurement, but 

(2) (2) 

leave it unspecified as M2 = {M^/ (g) N^,'} we then 

require a generalization of ([5]). Using the time-reversal 

(2) 

invariancc of the unitary interaction we have tr [M^7 «) 

N^VUaABU''] = ti[{eaABQ'')UeMlV <»N^^Ve'<U^, and 
from this we see that, for the pairing trick to work, the 
POVM elements of A^2 must themselves be valid quantum 
states of the same form prepared by 1 and the set of 
elements should be closed under the time-reversal oper- 
ator O. This on its own is a highly restrictive condition, 
and explains why forming a theorem for more general 
POVMs than the projective case is difficult. 



Application to semi-classical maximum work theorem. 
The above results, and in particular (|14p find simple ap- 
plication in a semi-classical maximum work theorem sce- 
nario [isj in which a quantity of ordered energy is ex- 
tracted from a primary quantum subsystem A. The pri- 
mary system is free to dump entropy in the form of heat 
into a heat sink B, with fast relaxation times, and ex- 
change mechanical work with a third (classical) adiabatic 
system C. 

On the assumption of conservation of energy for the 
composite system ABC and the adiabaticity of C the 
averaged relation leads to 

dWc <-dUA + TgdSA - Tb (dl) (15) 

where TAdSA '-^ (dq) corresponds heat flowing into A, 
and we assume for simplicity that no net work is done 
on A. The main point of this application is to display 
the contribution that the initial correlations between the 
primary subsystem and the reversible heat sink provide 
to the usual maximum work theorem, and in the process 
illustrate the well-known energetic value of correlations 

[IMS- 

Conclusion. We have extended the Jarzynski-Wojcik 
Exchange Fluctuation Theorem to the situation where 
we drop the assumption of molecular chaos, and allow 
correlations to exist in the composite state. These corre- 
lations results in a modification of the XFT relation and 
can enhance the probability of heat flowing in the back- 
ward direction. We have applied our results to deriving 
a semi-classical maximum work theorem for correlated 
systems. Our work highlights the difficulty of obtain- 
ing further results for situations without initial and final 
measurements of energy. 
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